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1 INTRODUCTION 

Since the first strong; gravitational lens was discovered in 
1979 (1 Walsh. Carswell fc Wevmannl) . hundreds have been 
found and studied (|Tre ul l20ld ). As well as serendipitous 
discoveries, dedicated surveys have sea rched for lens candi - 
dates. The Cosmic Lens All-Sky Survey |Browne et alj20o3 l 
identified 22 new strong gravitational lenses from imag- 
ing in the radio waveband. The Slo an Lens ACS Survey 
jBolton et al.ll200d ; lAuger et alJl200Sft has so far confirmed 
85 new gravitational lenses using spectroscopic data from 
the SDSS to identify candidates for follow-up with high- 
resol ution HST imaging. The Cosmol ogical Evolution Sur- 
vey (|Faure et all |200S| ; IJacksonl 120081 ) identified 67 lenses 
from visual inspection of HST i maging. In addition there 
are pr ojects, such as MU SCLES (|jackson et al.ll201ll ) and 
SQLS (inadaet alj|2012h . that are looking for new lenses by 
reprocessing data from large surveys such as UKIDSS and 
SDSS. 

There has been much interest in studying gravitational 
lenses because they magnify sources that would otherwise 
be too faint or distant to see. Historically these lens sys- 
tems have been difficult to identify in sufficiently large num- 
bers to be statistically useful. Some lens searches involve 
painstaking 'by eye' analysis of each object observed. Now 
new methodologies are being developed to improve the ef- 
ficiency of finding lenses and to improve the conversion 
rate from candidate to confirmed lensed source. These in- 
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The identification of strong gravitational lenses in large surveys has historically been 
a rather time consuming exercise. Early data from the Herschel Astrophysical Tera- 
hertz Large Area Survey (Herschel-ATLAS) demonstrate that lenses can be identified 
efficiently at submillimetre wavelengths using a simple flux criteria. Motivated by that 
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systems which have been, and will be, found by the Herschel- ATLAS. Analytical mod- 
els of lens statistics are tested with the current best estimates for the various model 
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lens density profile, for which we consider the singular isothermal sphere (SIS) and 
the Navarro, Frenk & White (NFW) approximations. The five lenses identified in the 
Herschel-ATLAS Science Demonstration Phase suggest a SIS density profile is pre- 
ferred, but cannot yet constrain Q\ to an accuracy comparable with other methods. 
The complete Herschel-ATLAS data set should be sufficient to yield competitive con- 
straints on Ha ■ Whilst this huge number of lenses has great potential for constraining 
cosmological parameters, they will be most powerful in constraining uncertainty in 
astrophysical processes. Further investigation is needed to fully exploit this unprece- 
dented data set. 

elude using colours dOfek et al l [20071 ), spectroscopic red- 
shifts || Auger et all 120091) , and other imaging processing 
techniques ((Marshall 200J). 

For a rigorous statistical study of gravitational lenses, 
they should all come from a single survey. This is because 
the probability of identifying lensed sources varies between 
surveys due to the different instrumental sensitivities and 
selection effects. 

For sources observed in the submillimetre waveband it 
has long been thought that galaxy number counts should 
fall off rapidly with increasing flux density. This is because 
of the rapid evo lution of proto-spher oidal galaxies into star- 
burst galaxies (|Granato et al.|[2004l ). This drop in counts 
with flux using submillimetre observations was anticipated 
to yield an effective way of identifying lenses (|Blainl Il996l ; 
iNegrello et alj|2007h . It has since been succe ssfully demon- 
strated using the Herschel- ATL AS data (|Negrello et al.l 
[2O10l ; iGonzalez-Nuevo et al.ll2012h . Lensed sources are iden- 
tifiable much more easily because the lensing magnification 
pushes them into higher flux densities which have intrinsi- 
cally low source populations. 

The Herschel Space Observatory, launched in 2009, is 
the only space-based observatory to cover a spectral range 
from the far infrared to the submillimetre. It provides a 
unique window through which to study large scale struc- 
ture. Th^^ersc/ieZAstrophysical Terahertz Large Area Sur- 
vey (|Eales et alj|20ld Herschel-ATLAS) will cover 550 deg 2 
of the sky making it the Herschel survey with the largest 
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area. iNegrello et all ||201Ch showed that strong gravitational 
lenses could be identified with almost 100% efficiency by se- 
lecting s ources which have 500/j n i flux densities greater than 
lOOmJy. iGonzalez-Nuevo et al.l (|2012l ) further improved on 
this technique by refining the flux selection to incorporate 
other Hersche l wavebands and optical/infrared data from 
the VIKING l|Fleuren et al.l 12013 ; [Sutherland et al.l l2012h 
survey to look for counterparts. This techn ique is forecast 
to yield ten times more lenses than that by INegrello et all . 
The expected e fficiency of 70% is lower that that achieved by 
INegrello et all but still impressive. Using these techniques, 
the Herschel- ATLAS team expect to find as many as a thou- 
sand strongly lensed candidates. 

We emphasise that this article considers the likelihood 
of strong gravitational lenses which are defined as lensed sys- 
tems with multiply imaged sources. Recent art icles on the 
effect of lensing on the so urce number counts dLima et al.l 
|2010| ; IWardlow et al.|[2012l ^ do not consider strong lensing in 
the sense that sources are multiply imaged, just that they 
are highly magnified. Note, this results in the lens probabil- 
ity distribution assuming a NFW density profile becoming 
much more similar to that of the SIS density profile than is 
seen in our analysis. Looking to the future, 'strongly' lensed 
sources are being selected by their excess submillimetre flux, 
or effectively by their magnification, with limited optical 
follow-up to confirm multiple imaging. Since requiring opti- 
cal follow-up may limit the number of confirmed sources, it 
may in future turn out to be statistically favourable to relax 
the multiple image requirement to gain a statistical advan- 
tage. However our approach, looking at strong gravitational 
lens systems with multiple images, is used t o be consistent 
with t he nature of the sources identified by INegrello et al.l 
(|2010h . 

The statistical properties of gravitationally 
lensed sources were fi rs t co nsidered in the 1980s 
l|Turner. Ostriker fc Gotd 1 1984 ) and later applied 
as a possib le way of c onstra i ning the cosmologi- 
cal model llLi fc Qstrikerl |2002|; iMitchell et al.l 120051 ; 
IZhang. Cheng fc Wul I2009T ; iDobke et al.l I2009T ). With the 
prospect of a thousand lenses being identified by Herschel- 
ATLAS, we take a fresh look at the analytical theory behind 
predicting strong gravitational lens statistic s in Sections [2] 
In Section [3] the lenses identified by INegrello et al 



are compared against the analytical predictions, and 
then in Section \5\ we consider the parameter constraints 
possible with the full Herschel- ATLAS data set. 

Unless otherwise stated we use cosmological parame- 
ters: n m = 0.2 7 Oa = 0.73, n = 046, h = 0.7, n = 0.97, 
and <7 8 = 0.81 (|Komatsu et alJboill) . 



2 LENS STATISTICS: THE THEORY 

A source is defined to be strongly gravitationally lensed 
when it has been multiply imaged due to the gravitational 
effects of mass between the source and the observer. The 
probability of a source being lensed is known as the optical 
depth (r). It depends on three main components: the cos- 
mology which determines the co-moving volume element at 
a given redshift; the normalised mass function (n) which de- 
scribes the number density of halos; and the lensing cross- 
section (a), which is the area in the lens plane where the 



separation between the lens and source is sufficiently small 
for strong lensing to occur: 



dz. 



dD, , \ 3 / \ / \ 

- (1 + 2[J n(m, ZjJ aym, z^) dm, (1) 



where D l is the angular diameter distance to the lens, z x is 
the lens redshift, z is the source redshift, and m is the virial 
mass of the lens (Schneider. Ehlers fc Falcolll993r i. This is 
explained in detail in the Appendix, where we describe how 
to calculate each of the different components. As well as 
the dependence on the lensing redshift z v the dependence 
on other parameters such as the total magnification of the 
lensed source or the angular separation of the images can be 
calculated. 

When comparing the theoretical predication for the dis- 
tribution of strong lenses with observational surveys, an- 
other factor that needs to be taken into ac count is magnifi- 
cation bias l|Turner. Ostriker fc Gottlll984 ). 



2.1 Magnification Bias 

Magnification bias is a correction to account for that fact 
that lensed sources are magnified and are consequently ob- 
served to be brighter than the population from which they 
are drawn. Therefore sources which are lensed come from 
a population with a lower flux. Since there are generally 
more fainter sources than there are brighter sources, mag- 
nification bias almost always increases the chance of finding 
a lens compared to if the lenses came from the same flux 
population at which they are observed. 

The magnification bias at a particular flux is calculated 
as the ratio of the number of lenses which will be magnified 
to that flux over the number of objects which intrinsically 
have that flux times the probability of lensing: 



B(F) 
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where P(fJ,) = t(/j,)/t is the probability of a magnification 
given that the source has been lenseqj. 

The overall bias is then calculated as the average of Eq. 
[2] over all fluxes: 



B = 



F ,. 1 x ' 
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/ N(F)dF 
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where F n is the flux limit of the observational survey. 

Generally the flux distribution is assumed to take the 
form of a simple power law N(F) = No (F/Fo)' 13 and so 
Eq. [3] then reduces to B — J r(/x)) dpi,. However, for this 
particular survey the slope f3 is so steep that pP grows 
faster than r(/i) decreases and so this integral does not 
converge. But assuming that this power law approximation 



1 Note that we do not he re incl ude a £t 1 factor 
as in iTurner. Ostriker fc Got i] jl984) . A s detailed in 
iFukugita fc Turned ||2007|) since both r and t(/z) are con- 
sidered in the same plane this 'correction' factor is not required 
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for N(F) holds for infinitely small fluxes is clearly incor- 
rect; there will come a point where the flux counts fall off. 
To account for this we introduce an exponential term that 
stops the growth of N(F) below a given flux F c meaning 
that N(F) will return to zero as F tends to zero. We use 
N(F) = N (F/F )~ fi exp(-(F/F c ) 2 ) and the magnification 
bias becomes: 



B = 



/ P( M ) / / F~P ex P (-(F/ M F c ) 2 ) dFdv 

oo 

/ F~? exp(-(F/F c ) 2 )dF 



(4) 



In this work we simplify the above expression by setting J* 1 ,, 
to 0. Eq. Uthen reduces to: 



B = j P{n)l*dp, 



(•>) 



which intuitively makes sense as it is just the average 
magnification. This work assumes that F.. = 0. 

Whilst Eq. [T]is relatively simple, it is more complex to 
implement in practice due to the uncertainties associated 
with each of the constituent components. In the rest of this 
section we explore the current best estimates for each of the 
components and consider the likely impact of uncertainty in 
each of them on the probability distributions for observing 
gravitational lenses. 



3 LENS STATISTICS: DEMONSTRATING THE 
THEORY 

In this section the different components of the lensing optical 
depth are considered. We discuss the current best estimates 
for each component and demonstrate the corresponding ef- 
fects of uncertainty on the optical depth. 

3.1 Dark Energy Density 

The community is, in general, confide nt in the ACDM cos - 
mologi cal model refined fro m CMB l|Bennett et al. 2003), 
BAO (|Percival et all l20ld) and supernovae IjRiess et al.l 
1998) observations. However it is still worth considering 
how the different cosmological parameters would affect these 
lensing statistics. If variations in the underlying cosmology 
significantly affect the lensing optical depth then the statis- 
tics of gravitational lenses could provide another comple- 
mentary test of the cosmology. On the other hand, if the 
cosmological parameters have very little impact on the lens- 
ing statistics then we can be confident that any assumptions 
we draw are not dependent on the accuracy of the ACDM 
assumption. 

The dark energy density (£2a) contributes to the lens- 
ing optical depth (described in Eq. [T] and the Appendix) 
by increasing the comoving volume the light from a source 
travels through. This makes it more likely that a source will 
encounter an object and be lensed. If we assume a flat uni- 
verse, this increase in Q,a is balanced by a reduction in the 
matter density (il m ). This decrease in il m will reduce the 
density perturbation growth rate and hence decrease the 
mean mass in the mass function. Figure [T] shows an example 
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Figure 1. The differential lensing optical depth as a function of 
lens redshift for Q\ = 0.66 (purple dashed line) and 0.7 (blue solid 
line). Thi s assumes a source r edshift z s = 2, SIS lens density pro- 
file and a lTinker et al . (2008) mass function. The increase in the 
comoving volume element and decrease the density perturbation 
growth rate combine so that the overall optical depth decreases 
with increasing Qa- 



of the differential lensing optical depth as a function of lens 
redshift for two different dark energy densities (Qa = 0.66 
and 0.73). The increase in the comoving volume element and 
decrease in the density perturbation growth rate combine so 
that the overall optical depth decreases with increasing Qa- 
In this example a 10% decrease in Qa results in a factor of 
two increase in the differential lensing optical depth. There 
is also a slight bias to lower z, with increasing Qa- The shape 
of the differential lensing optical depth does not significantly 
change as a function of 9 and fi (not shown) . 

Previous work has used results from the CLASS 
|Browne et al.l [20031 survey to constrain Q A . The CLASS 
data set, which includes 13 lenses, was used to c onstrain the 
total dark energy density Oa for a flat Universe (|Chae et al.l 
120021 : iMitchell et al.ll2005l) and the dark ener gy equation of 
state parameter 1 Zhang. Cheng fc Wull2009l ). For illustra- 
tive purposes, here we consider just varying the dark energy 
density for a flat cosmology. There are many other cosmo- 
logical parameters which could be constrained such as the 
curvature. This would affect the apparent angular separa- 
tion of the lens images. A positive curvature would make 
the angular sep arations appear sm aller than they would in 
a flat universe IjPark fc GottJI 19971 ). With the large number 
of lenses anticipated in the future it may also be possi ble to 
constrain modified gravity (|Zhao. Li fc Kovamall20l"ll ). 



3.2 Lens Distributions 

Understanding the lens distribution is important as lensing 
is a gravitational effect and can tell us about the distribution 
of all matter not just baryonic matter. The distribution of 
matter over-densities, or halos, is frequently approximated 
using an analytical approximation called the mass function. 

The mass function n(m, z) (sometimes written as 
dn/ dm) predicts the differential number of halos of a given 
mass as a function of redshift. It is frequently approximated 
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using the approach bv lPress fc Schechterl (|l974h : 



m 2 n(m,z) dm „, . dv 

^ — - — = vf(y) —, 

p m v 



(6) 



where p is the average matter density, v(m,z^) = 8 c /a , 
8 C is the critical density required for spherical collapse and 
a 2 is the variance in the linear density fluctuation field 
(see also Appendix IA3|) . Th ere have been a number o f es- 
timates proposed for f(v). iPress fc Schechterl (|l974l ) de- 
rived an approximation based on theoretical arguments. 
This was later shown to under/over predict the abundance 
of halos of high/low m ass relative to N-body simulations 
llSheth fc Tormenlll999l). Improved estimates hav e been pro- 
200ll; jWarren et al.ll200fj ). most 



posed (Jenkins ct al. 
cently bv lTinker et al 



( 20081 ) whose fit introduced a specific 



dependence on redshift. 

So how does uncertainty in the mass function n(m, z t ) 
affect the differential lensing optical depth? Figure [5] 
shows the differential le nsing optical depth for both the 
IPress fc Schechterl l|l974) and ITinker et all J200ct) app roxi- 
mations for the mass function. The lPress fc Schechterl mass 
function esti mates lower values of dr/dz l than that the 
I Tinker et al.l mass function. The shape of dr /dz x tends 
to favour: larger z t for the Pres s fc Schechterl mass func- 

mass function; and 
mass function cona- 



tion compared to the ITinker et all 
smaller 9 for the 



Press fc Schechterl 
— 1 



pared to th e Tinker et al.l mass function. This is because the 
IPress fc Schechterl mass function predicts a lar ger n at small 
m, and a smaller n at high m, compared to the I Tinker et al.l 
mass function. Integratin g over all m result s in a smaller 
dr/dzj and dr/dO for the lPress fc Schechterl mass function 
because the higher mass halos contribute most to the lensing 
optical depth. dr(p)/dp maintains the same shape for both 
mass functions (not shown). 

Whilst it is theoretically possible to generate a lensed 
source with a lens made purely of dark matter, the way in 
which candidate lens sources are confirmed relies on follow- 
up imaging to identify a lens candidate. Without an ob- 
served lens candidate it is likely the source would be re- 
jected as a potential lens. This is not something that has 
been accounted for in this work. The lensing optical depth 
is approximated with a mass function for the total matter 
distribution. This is worth noting as a potential source of 



3.3 Halo Density Profile 

Eq.rjjis most sensitive to the lens density profile component. 
The simplest model for a halo density profile that can ex- 
plain the flat rotation curves of galaxies is a singular isother- 
mal sphere (SIS). This has a density profile given by: 



p(r) = 



27rGr 2 



(7) 



where v is the velocity dispersion (see Appendix I A5[) . This 
profile has been successful in reproducing a num ber of fea- 
tures of gravitational lensing l|Auger et al.ll2009h . However 
the profiles of dark matter halos generated via numeri- 
cal simulations are better fit by the profile proposed by 
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Figure 2. The differential lensing optical depth as a function 
of redshift for z g = 2 and a SIS lens profile. The green dashed 
line is for the iPress & Schcchtei] mass function, which is com- 
pared to the pTinkeT'et al.l mass function (blue solid line). The 
IPress fc Schechterl mass fun ction underes timates the lensing op- 
tical depth compared to the ITinker et all mass function. 



lNavarr~ Frenk fc White] (| 19971 . NFW): 
p(r) = 



fe)( 



1 + 



(8) 



where r s is the scale radius and p s is the density at r a (see 
Appendix I A6|) . Neither the SIS or NFW density profiles is 
probably realistic for a lens made of a dark matter halo 
and baryon component. Here we use them for illustrative 
purposes to demonstrate the scale of uncertainty and its 
significant effect in the resulting optical depth. 

SIS profiles produce more strong gravitational lenses 
than the NFW profile. This is because they have a steeper 
inner density gradient (p oc r~ 2 as opposed to p oc r^ 1 
for the NFW) and steeper inner profi les are more effectiv e 
at producing multiply imaged lenses |Li fc Ostrikerll200l l. 
NFW profiles produce lenses with higher mean magnifica- 
tions than those for the SIS . This is because lenses that are 
less concentrated (such as the NFW) require the source and 
lens to be very closely aligned for strong lensing to result. 
It is this close alignment that results in a higher magnifica- 
tion (|Wvithe. Turner fc S pcrgcl 200l]). In other words lenses 
that are more concentrated (such as the SIS) do not need the 
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lens and source to be so closely aligned for multiple imaging 
to occur and this results in many more lower magnifications. 

Figure [3] compares the differential optical depth for the 
SIS and NFW halo profiles. As well as demonstrating the 
known results for dr/dfi discussed above, other effects are 
seen for dr/dz l and dr/dO. The NFW profile tends to pre- 
dict more small-separation lenses and fewer large-separation 
lenses than the SIS profile. Also the mean lens redshift is 
much smaller for the NFW density profile than the SIS. 

These results confirm that the theoretical strong grav- 
itational lensing probabilities are highly sensitive to the 
density gradient ne ar the centre of dark matter halos 
l|Li fc Ostrikerl 120021). This is an area of tension in cur- 
rent observations l|Comerford et alJliooil ; IShu etTaJ feoOS). 



High resolution observations of galaxy rotation curves sug- 
gest halo density profiles with a flat central core. However, 
studies of gravitational lensing statistics sugg est a steep, 
or 'cu spy', density profile for halo cores l|Chen fc McGaughl 
|200S| ). 

A potential resolution for this is that the profiles 
of the gravitational lens populat ion are biased wit h re- 
spect to the general population ( M andelbauml [20091 ') ; al- 
ternative profiles have also been suggested. Following 
iNavarro. Frenk fc White] (Il997f), h igher resolution numeri- 
cal simulations ([Moore et al.l 19991 ') suggested that the cen- 
tral density gradient should lie somewhere in between that 
of the NFW (p oc r" 1 ) and SIS (p oc r~ 2 ) profiles. It may 
be that par ameterising the slope of the density profile with 
an Ein asto ( Retana-M ontenc gro fc Frutos-Alfaro 20 111 ) or a 
Sersic |Mandelbaumll2009l ; lEh'asdottir fe Mollerj 120071 ) pro- 
file would be more appropriate. 

Alternatively, CLASS found a number of lenses at low 
separation, but none at higher 8, which suggested a combina- 
tion of NFW and SIS density profiles. They proposed that if 
the mass of the lens is greater than a given cut-off m c then 
the lens density profile is approximat ed using a NFW ap - 
proximation, else a SIS profile is used (|Li fc Ostrikerll2002h . 
This was not motivated by direct observations/simulations 
of halo profiles, but it has been shown to provide a rea- 
sonable fit to lensi ng statistics for a mass cu t-off scale of 
m c = lO^Mo/i" 1 (|Zhang. Cheng fc Wdl2009T ). 

The NFW halo density profile could well be accurate 
but the halo concentration could be incorrect. The concen- 
tration of the halo (defined as c = J'VirA's) has been stud- 
ied extensively over the last decade. On average, the me- 
dian halo concentration decreases with increasing mass and 
redshi ft; this can be approximated us i ng analytical expres- 
sions (iNayarro, Frenk fc White! Il997l ; iBullock et all l200ll ; 
lEke et al.ll200lh . 

All these suggestions assume a spherically symmetric 
lens, but realistically the halos are more likely to be irreg- 
ularly shaped with a spherically symmetric average. Some 
way of accounting for this variation may well be required to 
provide a consistent explanation across dark matter simula- 
tions and lensing studies. 
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Figur e 3. The di fferential lensing optical depth for z s = 2.0 
and a iTinker et al.l mass function. When not varying the lensing 
redshift is fixed to = 0.5. The solid blue line shows the SIS 
compared to the NFW (red dashed line) lens profiles. The SIS 
profile produces more lenses than the NFW profile because it has 
a steeper inner density gradient. 



3.4 Source Distributions 

The analytical calculation of the lensing optical depth in Eq. 
[T] only describes the probability of a source being lensed. It 
does not account for the probability of a source being at a 
particular redshift in the first place. So if you wanted to esti- 



mate the total optical depth this source redshift distribution 
needs to be measured observationally. 

For the Herschel sources it is generally not possible to 
measure accurate spectroscopic redshifts due to the lack of 
emission lines in that band and the low angular resolution. 
Much effort has been put into deriving accurate source 
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ID 


h 






9 


67 g±0.057 


1 577±°- 008 




11 


.7932± 0012 


L786 ±0.005 




17 


0.9435 ±0 - 0009 


2.308± - 011 




81 


2999 ±o.ooo2 


3042 ±o.ooi 


25±7 


130 


0.2201± 002 


2.6260± 0003 


6±1 



Table 1. Estimates of the source (z s ) and lens (a.) redshifts and 
magnification (ii) for the five lens es discovered in the Herschel- 
ATLAS SDP jNeerello et al.l2010h . The quoted errors correspond 
to a 68% confidence interval. 



redshift distributions for the Herschel source catalogues via 
other methods such as matching so urces to optical cou nter- 
parts with spectroscopic red shifts dSmith et. all [201lh . CO 
line spectroscopy follow-up (|Harris et al.| 2012), and using 
models based on average galaxy SEDs l|Amblard et al.ll2010l; 
lLapi et al.l l201ll ; iPearson et"afl 120121 '). All these methods 
have their own biases/constraints, so a potential test that 
could be done with this lens work would be to fix all the 
other components and investigate the source distribution. 



In this section we have seen that the lensing statistics 
are sensitive to both astrophysical and cosmological param- 
eters. The lensing optical depth is most sensitive to the lens 
density profile, both in magnitude and as a function of all 
z,, n and 6. The mass function has little effect on the optical 
depth as a function of /x, but does for z, and 9. Of all the 
parameters considered, the optical depth was least sensitive 
to the dark energy density. It should be noted that before 
competitive constraints on the cosmological parameters can 
be achieved, the uncertainty in the properties of the lensing 
galaxies needs to be constrained. 



4 LENS STATISTICS: HERSCHEL- ATLAS SDP 

Five new lensed sources were discovered using the flux crite- 
ria technique on the H erschel-ATLAS Scien ce Demonstra- 
tion Phase (SDP) data l|Negrello et al.ll2010l ). Table [1] sum- 
marises the estimated source and lens redshifts for those 
lenses. Whilst this is probably not enough lenses with which 
to do convincing statistics, we can still consider whether 
or not these observations are consistent with the analytical 
models discussed in Section [2] 

Here we consider both the absolute differential optical 
depth as a function of lensing redshift drizA /dz^ and the 
conditional differential optical deptrfl 



1 dT(u) / \ 1 dr(z ) 



which describe the probability of a lens being at a certain 
/i or z l given that lensing has occurred. We do not consider 
P(0) here because estimates are not yet available for the 
image separations of the Herschel- ATLAS SDP lenses. 

2 Note here we consider drifi)/dfi = [dr(> [i-\-dfi) — dr{> fj,)]/dfj, 
not g(t(> /J.)/dfi as in Section [3] 



By choosing to consider the conditional probability, 
rather than the absolute probability, the calculation is sim- 
plified. The magnification bias (see Section 12. ip does not 
need to be calculated as the conditional probability is effec- 
tively the 'normalised' optical depth for fixed z a . The con- 
ditional probability is therefore the expected distribution of 
the lenses as a function of fi or z, given that lensing has 
occurred. 

We start by looking to constrain the lens density profile 
since this is the component to which the statistics of strong 
lenses appeared to be most sensitive in Section [231 We start 
by considering the conditional probability as a function of 
magnification P(/i) and redshift P fz,) to constrain the lens 
density profile between the NFW and SIS. 



4.1 Constraining the Lens Density Profile with 

Magnifications are availa ble for the Herschel- ATLAS SDP 
lenses ID081 and ID130 (|Negrello et alj 1201(f ). Theoretical 
distributions of the conditional probability as a function of 
magnification were calculated as d escribed in the Appendix. 
They were calculated assuming a iTinker et al.l mass func- 
tion, and a standard ACDM cosmology. The source and lens 
redshifts are as per the estimates in Table [1] The resulting 
P(/x) are shown in FigureUfor the two different estimates of 
the lens density profile (SIS and NFW), along with the two 
observational estimates of the magnifications for ID081 and 
ID130. This shows that the magnification of ID081 source is 
within the 99% confidence internal of the theoretical prob- 
ability found for both the SIS and NFW density profiles. 
The magnification for the ID 130 source is consistent with 
the model with the SIS density profile at a 95% level, but 
inconsistent with the model with the NFW profile at over 
a 99% level. These results appear to favour the SIS density 
profile. 

One way of quantifying the combination of the ID081 
and ID 130 results to choose between the NFW and SIS op- 
tions is to use a likelihood. The likelihood of the data given 
the model parameters is defined as: 



£ (^l^ S i S/N pw) =n p (^)' 



(9) 



where N, = 2 is the number of lensed sources being consid- 
ered. To give the most favourable result for each of the lens 
models, the highest probability associated with the range of 
observed magnification is chosen for each P(/ii). Calculat- 
ing this gives £(/ii | p sxs ) = 3 x 10 3 and C[fii | P NFW ) < 
5x 10~ 6 which confirms that a SIS density profile is preferred 
overall. 

The observational estimates for the magnification in 
Figure [3] are consistently higher than those preferred by the 
SIS lens density profile. We note that magnifications close 
to 1 will be much harder to identify observationally than 
higher magnifications. Therefore we would expect that the 
observational values would be biased against magnifications 
close to n — 1. This suggests that the SIS likelihood is ac- 
tually higher. Overall we see the Herschel- ATLAS SDP lens 
magnifications show that the SIS density profile is much pre- 
ferred over the NFW density profile. 
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Figure 4. Conditional probability of finding a lens at a certain 
magnification given lensing has occurred, for source and lens red- 
shifts in Table [T] Solid blue lines use the SIS lens density pro- 
file; dashed red lines use the NFW lens density profile. Shading 
indicates 68% confidence intervals. Solid green vertical lines in- 
dicate the estimated magnification of gra vitational lenses ID081 
and ID130 found bv lNegrello et alj fcOloT) . 



4.2 Constraining the Lens Density Profile with 

In Section [3] we saw that the lens density profile strongly 
impacts P (z.) and P(0) as well as P(fJ-)- Data is not yet 
available for the angular separation of the lensed images, 
therefore we just look at P (zj) here. There were only two 
data points available for the lens magnifications but lens red- 
shift estimates are available for all of the Herschel-ATLAS 
SDP lenses. 

Figure [5] shows P(z l ) for all five of the Herschel-ATLAS 
SDP lenses. Again P(z l ) was calcula ted as described in 
the Appendix assuming a lTinker et alj mass function and a 
standard ACDM cosmology, for any angular separation and 
magnification (i.e. 9 > and /j, > 1). The source redshifts 
are as per the estimates in Table [1] 

Three of the five Herschel-ATLAS lensed sources 
(ID009, ID011, and ID017) show lens redshifts most con- 
sistent with 68% confidence limits for the SIS density profile 
model. These positions are still consistent with the NFW 
model at around the 95% confidence level. The other two 
Herschel-ATLAS lensed sources (ID081 and ID130) are both 



at much smaller lens redshifts. The z l are consistent with 
both the SIS and NFW models again at a 95% confidence 
level. 

To quantify these results we can again calculate likeli- 
hood as: 



1 % ^SIS/NFW 



Calculating this gives Cfz^ | A™) 



(10) 



0.014 and 



0.026. 



So, in contrast to what was 
shown with P(fJ,), the NFW profile is preferred over the SIS 
density profile. 

Again we should note that there is likely to be an ob- 
servational bias here. The way in which these lens systems 
are selected requires a lens object to be identified observa- 
tionally. Objects which are closer will tend to be brighter 
and therefore easier to identify. Therefore it is likely that 
there is a selection bias for nearer lenses which has not been 
taken into account. This bias favours the NFW profile, so if 
it were taken into account the SIS profile would likely have 
a more favourable likelihood. 



4.3 Constraining the Dark Energy Density with 

dr fdz l 

Here we consider whether the dark energy density can be 
competitively constrained using the lens redshifts z } of the 
Herschel-ATLAS SDP lenses as compared to the theoretical 
dr^Zj) /dz r The tests above provide conflicting evidence for 
the SIS and NFW lens density profiles. Later in Section [5] 
we will put these results into context using simulations, and 
project how many lensed sources are needed for a definitive 
result. Meanwhile, for illustrative purposes, here we assume 
a SIS density profile. Similarly since there is no information 
on the angular separation, we make the ass umption that the 
mass function is well fit by a iTinker et alj mass function. 

Now we consider absolute probabilities we also need to 
take into account the other sources that were not lensed. 
This is to account for the probability of lensing occurring as 
well as the distribution of the lenses as a function of z l . The 
likelihood is given by: 



i t \ ™ 



r («si)), 



(11) 



2=1 



where N u is the number of unlensed sources. Given an es- 
timated redshift distribution Af(z s ) then, because the op- 
tical depth is small, we can approximate C^z^ j QaJ as 
l|Mitchell et al.ll2005T ): 



In C [zji j Qa) = E In 

i=l 

if. 



i=l 



dr 




Z x i) 




dz. 




(It 




hi) 




dz. 




dr 




Zyi) 


dz, 



£>(l-r(ay)) 

2=1 
2=1 



N(z a )T(z a )dz s 



z B i and Zj< are the source and lens redshifts of the individual 



8 Short, Pearson, Coles & Eales 



ID009 ID011 




Figure 5. Conditional probability of finding a lens at a certain lens redshift given lensing has occurred, for all five Herschel-ATLAS 
lenses. Source and lens redshifts are in Table [T] Solid blue lines use the SIS lens density profile; dashed red lines use the NFW density 
profile. Shading indicates 68% confid ence intervals. Solid green vertical lines indicate the estimated lens redshifts of the gravitational 
lenses found bv lNegrello et aL 



Herschel- ATLAS SDP lenses but z s in the second term is a 
discrete range from the minimum to maximum redshifts of 
Herschel- ATLAS sources observed. 

The redshift distribution of the Herschel-ATLAS 
sources is quite uncertain; for illustrative purposes here we 
assume that all the sources are at a redshi ft of 2. We take 
the number of unlensed sources as 252 from lClements et al.l 



(2010.); the resulting Q,a is invariant to small changes in 
7V U . Figure [6] shows the resulting likelihood as a function 
of the dark energy density (assuming a flat universe). This 
corresponds to an estimate of Qa = 0.81 j^'j^. This is 
consistent with the current best estimate for £7a at the 95% 
level. 
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Figure 6. Likelihood of the observed lens redshifts of the Her- 
schel- ATLAS SDP lenses given a range of dark energy densities 
C2a- The shaded regions represent the 68% confidence region. The 
vertical green line represents the current most favoured SIa- 



To summarise, there are some conflicting results, which 
is not surprising given the current limited number of 
confirmed gravitational lensed sources. In the next section 
we put these results in context using simulations. We also 
consider what constraints can be achieved with a larger 
number of lenses. 



5 LENS STATISTICS: FULL 
HERSCHEL-ATLAS SURVEY 

So far just five strong gravitational lenses have been con- 
firmed in the Herschel- ATLAS SDP data. With the full 
data set it should be possible to extract enough lensed 
sou rces with which to do more meaningful statistical analy- 
ses (|Negrello et alJl201CtrGonzalez-Nuevo et alj|2012h . Here 
we consider whether the full set of 100-1000 lenses will be 
sufficient to constrain the different models for the density 
profile, mass function and Qa discussed above. 

The methodology is described in more detail below, but 
here is a general overview. The differential optical depth was 
calculated as described in the Appendix using different den- 
sity profiles, mass functions and J1a- A thousand sets of N, 
lenses were generated by randomly sampling the differential 
optical depth at different redshifts z lt magnifications fj,, and 
angular separations 8. The corresponding likelihoods were 
then calculated as described in Equations [9] [10] and [11] for 
each of the sets. The set size N l is varied to see how success- 
fully the different approximations for p, n, and Qa can be 
distinguished. Unless otherwise stated the source redshifts 
are all simply assumed to be z s — 2. 

This work could be expanded to incorporate a full 
Bayesian analysis. It is not considered here as we do not 
consider a sophisticated enough halo model, and there is 
not yet enough data, to justify such an approach. 

First we look at constraining the lens density profile 
using P(fJ,)- Ignoring the other possibilities (as discussed in 
Section [3} here we focus on the SIS and NFW profiles. A 
thousand sets of = 2 lens magnifications were generated 
by randomly sampling P(n) for both SIS and NFW den- 



Figure 7. Normalised distribution of number counts of likeli- 
hoods realised from random sampling of theoretical probability 
distributions P(fJ.) for SIS (solid blue line) and NFW (dashed red 
line) density profiles. The plot is for sets of N, = 2 lensed sources; 
the top plot is for lens redshifts generated assuming a NFW pro- 
file, and the bottom plot for SIS. The shading represents 68% 
and 95% confidence intervals. The two distributions only slightly 
overlap so it is fairly easy to distinguish between SIS and NFW 
using this approach - it will be even easier with more results. 



sity profiles. The iTinker et alJ mass function and standard 
ACDM parameters were assumed. Two versions of the 1000 
sets were generated for each of the SIS and NFW density 
profiles. The associated likelihoods were calculated for each 
of these sets for both the SIS and NFW lens profiles as de- 
scribed in Eq. [9] So in all four different combinations were 
considered. For each combination, the number counts for a 
range of different likelihoods were calculated. These results 
are shown in Figure [7] 

The blue and red curves are for likelihoods calculated 
with SIS and NFW density profiles respectively. The top 
plot is for sets of lens sampled from P(fi) using a NFW 
density profile; the bottom plot uses P(p) for SIS density 
profiles. The two distributions only slightly overlap in the 
first instance, and not at all in the second case (the NFW 
distribution is at smaller L than on the plot). It is easy to 
distinguish between SIS and NFW density profiles because 
their likelihood are distinct even with this small amount of 
data. 

The actual likelihoods found from the SDP data 
above in Section [3] were log£(/Xi | p SIS ) = —2.5 and 
log £(/ii | P NFW ) < —5.3. These are both consistent with 
the bottom plot which gives likelihoods for lens redshifts 
sampled from P(fJ.) assuming a SIS density profile. 

Whilst in Section[4]P(/ii) was found to be more consis- 
tent with a SIS model than a NFW model, f (z^) suggested 
the opposite. Here the same test as above was performed 
again but this time sampling the lens redshift as opposed 
to the magnification. Figure [8] shows the results. The top 
plot is for sets of N, = 5 lensed sources (as we have for the 
current SDP sample). The top half shows the distribution 
of likelihoods found using SIS and NFW probabilities for 
lenses randomly sampled from the NFW P (/i). Similarly for 
the bottom plot but here the lenses are randomly sampled 
from the SIS P{fJ,)- The two distributions overlap so it is 
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Figure 8. Normalised likelihoods realised from ran dom sampling 
of theoretical probability distributions f (zj) given a lTinker et al.l 
mass function and standard ACDM cosmological parameters, for 
SIS (blue solid line) and NFW (red dashed line) density profiles. 
The top plot is for sets of N, =5 lensed sources (as we have for 
the current SDP sample); the lower plot is for N, = 20 lenses. 
The top half shows the distribution of likelihoods found using 
SIS and NFW probabilities for lenses randomly sampled from the 
NFW P(fJ-)- Similarly for the bottom plot but here the lenses are 
randomly sampled from the SIS P(fJ.)- The shading represents 
68% and 95% confidence intervals. The two distributions overlap 
for N, = 5 so it is less easy to distinguish between SIS and NFW 
profiles; but at N, = 20 they are fairly easy to distinguish. 



less easy to distinguish between SIS and NFW using P{z 1 ) 
than when using -P(/i). Likelihoods found above in Section[4] 
gave logZI^i | pgjg) = -1.9 and log^Zji | p NFW ) = -1-6. 
Although the NFW likelihood is higher we see from Figure[8] 
that is not inconsistent with a SIS lens. More data is needed 
to properly distinguish between them; N, = 20 lenses would 
be sufficient as seen from the bottom plot. 

Therefore to summarise, whilst the results for the den- 
sity profile appeared conflicted initially in Section^ here we 
see that actually they are both most consistent with a SIS 
density profile. 

Values of the angular separation of lensed images are 
not yet available, but here we consider whether using 9 
would allow us to constrain the best fit mass function. Fol- 
lowing the same approach as above but now sampling 9 from 
P(9), the results for N. = 100 lenses are shown in Figure [9] 
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Figure 9. Normalised likelihoods realised from random sampling 
of theoretical probability distributions P(6) given a SIS density 
profile and standard ACDM cosm ological param eters, for Press- 
Schecter (green dashed line) and iTinker et al.l (blue solid line) 
mass functions. The plot is for sets of JVj = 100 lensed sources. 
The top half sho ws the distribution of likelihoods found using 
I Tinker et al.l and iPress &: Schechterl m ass function probabilities 
for lenses randomly sampled from the ITinker et ail P{8). Simi- 
larly for the b ottom plot but here the lenses are randomly sam- 
pled from the IPress fc Schechterl P(9). The shading represents 
68% and 95% confidence intervals. The distributions are not dis- 
tinct so it is not possible to distinguish between the two mass 
functions using just the angular separations of 100 lensed images. 



This shows that even with angular separations for 100 lensed 
sources it would n ot b e possible to d istinguish between the 
IPress fc Schechterl and ITinker et all mass functions. Of the 
order 500 sources would be required as a minimum. The 
number of detections is expected to be this high, but the 
shear amount of follow-up that would be required to get 
enough angular separations 9 makes this approach unfeasi- 
ble. 

Finally, we consider constraining the dark energy den- 
sity with the absolute differential optical depth dr/dz 1 . 
The lens redshifts are randomly sampled from the optical 
depth assuming the standard ACDM cosmological parame- 
ters. Rather than considering two alternatives of Ha, here we 
consider a range of Qa to calculate the likelihoods. There- 
fore rather than comparing the normalised distributions of 
C(dr I dz.) , instead we take the m ean value C. f or each £Ia 
and this is plot in Figure [6] The ITinker et alj mass func- 
tion and SIS density profile are used. The plot is for sets of 
N x = 5 (as we have for the current SDP sample), 100 and 
1000 lensed sources. The shaded areas are 68% limits. 

We see that the SDP estimate of Qa is consistent with 
the sort of errors expected for a sample of 5 lens redshifts, 
but with ~ 100 lens redshifts the error on Qa should be of 
order 0.01 i.e. comparable to the current best methods. 



6 DISCUSSION AND CONCLUSIONS 

Until recently the number of strong gravitational lenses that 
had been identified was relatively small for any particular 
survey. Now this situation is about to change dramatically 
as future surveys are expected to yield numbers of strong 
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Figure 10. Mean likelihood as a function of dark energy density 
realised from random sampling of theoretical proba bility d istri- 
butions dr / dz^ given a SIS density profile and iTinker et al.l mass 
function. The shading represents 68% confidence intervals. As- 
suming the other astrophysical uncertainties can be constrained, 
we see that competitive constraints on Oa can be achieved with 
the full Herschel-ATLAS data set. 



lenses into the hundreds and thousands. One of these sur- 
veys, Herschel- ATLAS, expects to yield up to 1000 lenses. 
In this article we considered the statistics that can be done 
with this new data set of strong gravitational lenses. 

First we reviewed the current analytical theory. We 
noted that the standard method for calculating magnifica- 
tion bias was not appropriate here. This is because the fast 
fall in number counts with flux in the submillimetre wave- 
bands means that the standard calculation gives an infinitely 
large magnification bias. We therefore proposed a more real- 
istic method. This effectively means the magnification bias 
is equal to the mean expected magnification of the gravita- 
tional lenses. 

We saw that lensing statistics are sensitive to astro- 
physical properties, as well as the cosmological parameters, 
most significantly the lens density profile. Although so far 
only five lenses have been confirmed in the Herschel- ATLAS 
SDP data, limiting its the statistical value, we have shown 
that these can be used to constrain the lens density pro- 
file. Initially the likelihoods calculated from each of P(n) 
and P\Z\) suggested conflicting results between the SIS and 
NFW lens density profiles. But simulations of the expected 
uncertainty of these likelihoods due to the small number of 
lenses showed that the results are perfectly consistent with 
a SIS, not a NFW, density profile. Of course neither of the 
SIS or NFW density profiles is probably realistic to explain 
the complex structure of the dark matter halo and baryon 
component. Here they are used to illustrate the effect of the 
current uncertainty in the density profile on the resulting 
optical depth. 

We considered just one of the cosmological parameters 
for illustrative purposes. Whilst five lenses is not yet enough 
information to constrain the dark energy density, the full 
Herschel- ATLAS data set should be more than sufficient to 
provide constraints competitive with other methods. How 
the mass function could be independently constrained was 
more uncertain and further work is clearly required to un- 



derstand how this and other other uncertainties can most 
effectively be constrained. 

Here we have seen that the dominant uncertainties in 
gravitational lens statistics are not cosmological but astro- 
physical i.e. understanding the nature of the lensing objects. 
Whilst more investigation is required, we can see that new 
data sets of strong gravitational lenses will be fundamental 
for constraining both halo astrophysics as well as cosmolog- 
ical parameters. 
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APPENDIX A: THEORETICAL STATISTICS 
OF GRAVITATIONAL LENSES 

Al General Formulation 

A source is denned to be strongly gravitationally lensed when 
it has been multiply imaged due to the gravitational effects 
of mass between the source and the observer. The probabil- 
ity of a source being lensed is described by the optical depth 
(t). It depends on three main components: the cosmology 
which determines the co-moving volume element at a given 
redshift; the normalised mass function (n) which describes 
the number density of halos; and the lensing cross-section 
(a), which is the area in the lens plane where the separation 
between the lens and source is sufficiently small for strong 
lensing to occur: 

f f dD l 3 

t(z s ) = / / —j— -(1 + Zj) n(m, Zj) o~[m, zj dmdz v (Al) 
o o 1 

where m is the mass of the lens, z, is the redshift of 
the lens, z s is the redshift of the source, and D l is the 
angular diameter distance from the observer to the lens 
l|Schneider. Ehlers fc Faicol Il993l) . Below we describe how 
each of these components can be calculated. 



#1 




Figure Al. Illustrative diagram of a gravitational lens system. 
The source, lens and observer are positioned at points S, L and 
O respectively. The source and lens are at distances D s and D l 
from the observer, and separated by a distance D. . f and n are 
distances in the lens and sources plans from the observer-lens axis 
(dot-dash blue line). The angular separation of the source from 
the observer-lens axis, as would be seen by the observer without 
lensing, is 0. However the observer does not see the source directly 
but a modified image (I) of the source due to the deflection by 
the lens. This image appears to be at an angular separation from 
the observer-lens axis of 8. The actual deflection angle a, and the 
reduced deflection angle a, can be estimated using the thin lens 
approximation as detailed in Eq. IA9I 



A2 Angular Diameter Distance 

The angular diameter distanc43 between the two redshifts 
z a to Zb is given by: 



D ah = 



c 1 



E(z)dz, 



(A2) 



Ho 1 + z b 

z a 

where E(z) = [« m (z = 0)(1 + zf + Q A (z = 0)] ' 5 . 
A3 Mass Function 

The differential number density of halos (n, sometimes writ- 
ten as dn/ dm) is commonly e s timat ed using this expression 
derived bv lPress fc Schechterl d 19741 ): 



m 2 n(m, z ) dm 



dv 



(A3) 



p = ri m p c is the average matter density, v(m,z^) = 
5 2 /a 2 (m, »,), 5 C — 1.686 is the critical density required for 
spherical collapse and a 2 is the variance in the linear density 
fluctuation fielcjfl. Two ap prox imations for f(y) ar e those by 
IPress fc Schechterl |l974T ) and lTinker et all feoOSt ): 



/psM-Vi exp (^)' 

and 



f T (y) = A 



where A 



exp 



—vc 



(A4) 



(A5) 



0.186(1 + z)-°' 14 , a 



1.47(1 + z) 



2.57(1 + z)~ a , a = 1.068, and c = 1.19. 



3 This definit ion assumes a flat un iverse. 

4 We use the lBardeen et al] l|1986l ) transfer function to calculate 
the linear power spectrum. 



A4 Lensing Cross-section 

The gravitational lensing cross-section (a) is the area in 
the lens plane where the lens and source are within the 
critical angular separation I3 r f or strong lensing to occur 
^Schneider. Ehlers fc Faicoll 19931 ): 

a = n/3 2 D 2 . (A6) 

The critical angular separation /3 C that produces a depends 
on the shape of the lens density profile. Two different ap- 
proximations for the lens density profile (SIS and NFW) are 
considered in Sections lA5l and [A6l First the general approach 
for spherically symmetric lens density profiles is reviewed. 

The critical angle (/3 C in Eq, IA6|) can be estimated using 
the lens equation: 



(A7) 



where a = aD ls /D s (see Figure lAl) . Note this assumes a 
flat space time so that the Euclidean relation separation = 
angle x distance holds. 

The lens equation can alternatively be written in terms 
of the radial distances, r)/D s — f/-D, — aD ls /D s . This can 
be simplified further using dimensionless distances y = 77/770 
and x = £/£o where r/o = £oD s /D 1 , to give the lens equation 
of the form: 



V 



D. D. 

t ls 1 



(A8) 



Using the thin lens approximation, which uses the fact 
that the lens is thin compared to the total path traveled by 
the light, the deflection angle a can be calculated as: 



a(x) 



4Gm(x) 



(A9) 



x£oc 2 

where m(x) = 4ir f v f p(r)dzdv and r 2 /r s 2 = z 2 + x 2 . 

The distance in the lens plane subtended by the critical 



14 Short, Pearson, Coles & Eales 
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Figure A2. Source position (r\) as a function of lens position 
(£) from the lens equation (Eq. \k&\ for SIS (green solid line) 
and NFW (blue dashed line) profiles. The diagonal red dotted 
line shows the case where there is no lensing. This example is 
for z s = 3, Zj = 0.5 and m = 10 15 Mq . For the NFW lenses, 
strong lensing (generation of multiple images) occurs for source 
positions with magnitudes less than the local extrema. For SIS 
lenses, there are no local extrema but here strong lensing occurs 
for ?7 < |jj(0)|. 



angle is simply y c — P c D s /rio. This is called the critical 
radius and is calculated as the value of y{x) where dy/dx = 
0. For values of y below y c , multiple images are formed. We 
will see below that two multiple images are formed in the 
case of the SIS density profile, and three in the case of the 
NFW density profile. 

A5 Singular Isothermal Sphere (SIS) density 
profile 



The SIS density profile is given by: 

„,2 

p(r) = 



(A10) 



km/s, (All) 



2nGr 2 ' 

The velocity dispersiorjf] can be approximated as 

v(m, Zl ) = 92.3[A° u , 5 r (^)ig( 2l ) iol3 ^ 1Ms ^ 

where A vir (z l ) = 18tt 2 + 60\n(z } ) - 11 - 32[0(z 1 ) - l] 2 and 
fi(z,) = « m (0)(l + z^/E^) 2 (Mitchell et alj|2005h . 

The deflection angle (Eq. IA9|l becomes a = Atvv 2 /c 2 
and the lens equation (Eq. IA8|l reduces to y = x ± 1 where: 

4™ 2 D ls D x 



e 



D. 



(A12) 



The critical value of y below which strong lensing (or 
multiple images) occurs is \y c \ = 1. We see from Figure 
A2I that two images are generated in this case. The lensing 
cross-section is given by: 

2 



fl 2 n 2 TTVcVoD^ 

a sis = ^Pc^i = - — 



D 2 



4tiVL> L>, 



c 2 D, 



(A13) 



5 Note, this velocity dispersion is from SDSS estimates for the 
CLASS survey. 



Other properties of the lenses, such as the magnification 
and angular separation, can also be calculated. The proba- 
bility that a lens will be magnified more than p is included 
in the cross-section as: 



a{> n) = aA(> ft), 



(A14) 



where A(> /i) describes the fractional area satisfying mag- 
nifications greater than fi: 



A{>^) = — 2 / n v Q \p(y) -t4 d v> 



(A15) 



where O is a step function. The magnification of an image 
is given by: 



ydy 

x dx 



(A16) 



l|Schneider. Ehlers fc Faicolll993D and since the solutions of 
the SIS lens equation are x± — y±l (|y| < 1), this becomes: 



M± {y) 



± 1. 



So the total magnification p(y) — /-t-(y) + p+(y) 
Substituting this into Eq. IA15I gives: 



(A17) 
= 2/y. 



A(>p) = — 
nyf 



- 2 [ *yQ[- 
Vc J y 



■u] dy = — for /i > 2. (A18) 



Therefore o" SIS (> p) = 4o" SIS /it -2 for /i > 2. 

Similarly the probability that a lens will have images 
with an angular separation greater than 6 is included in the 
cross-section as: 



,(>0) = a SIS eo?-< 



(A19) 



The separation between the two images can be seen from 
the lens equation to be Ax — 2. The angular separation 
becomes $ = £ Ax/D = 8ttv 2 D /c 2 D s . 



A6 iNavarro. Frenk fc White! (NFW) density 
profile 

INavarro. Frenk fc White] (|l997l ) proposed a halo density 
profile of the form: 

Ps 



p(r) 



7-U + ^ 



(A20) 



where r s is the scale radius and p a is the density at r s . 
The scale radius r B can be calculated as r s = r ■ Ic where 
the virial radius (r . ) is defined as the radius at which the 
mean density of the mass is 200 times that of the mat- 
ter density. The con centration can be approximated using 
jBullock et al.ll200ll ): 



c(m ' 2 ' ) = TTz(i 



1 + z x V 1.5 x 10 13 M G 



(A21) 



and p B can be calculated using, m = 47rp s r s 3 [ln(l + c) 
c/(l +c)]. The lens equation (Eq. IA8I) becomes: 



y = x g(x 

x 



(A22) 
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where fo = r s , fi s = 4p s r s /£ crit . The critical surface density 
is given by: 

E crit =c 2 DJ4nGD l D ls , (A23) 
and: 

X oo 

g(x) = juj(u 2 + v 2 y°- 5 (l + (u 2 + v 2 )°' 5 )-' 2 dvdu. (A24) 



It can be seen from Figure [A"2l that the critical value of 
y below which strong lensing (or multiple images) occurs is 
when dy/dx = 0. So the cross-section is given by: 

a NFW ( Zs ,^,m) =ttP c D 2 =n(y^jD 2 = ny 2 c r s 2 . (A25) 

As for the SIS profile, the dependence of a on other 
factors such as magnification and angular separation can 
be calculated. The lensing cross-section can be calculated in 
the same way as described in Equations lA14l and lA15l In the 
case of the NFW profile this does not reduce into a simple 
expression as is the case for the SIS. Instead o" NFW (> /U) 
needs to be calculated using the full expression detailed in 
Equations IA14I to IA16I This calculat ion is extens i ve an d a 
good approximation is available from lOguri et all ((2002). 

Again to calculate the dependence of the cross-section 
of the angular separation, o" NFW (> 9) ~ <j nfw O('& — 8). The 
separation between the two furthest images can be seen from 
the lens equation to be Aa; = 2x c where x c = /3 C D S . The 
angular separation is therefore $ = £ > o&.x/D l — 2r s x c /D l . 



